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We present a theory which describes the transport properties of the interfacial region with respect 
to heat and mass transfer. Postulating the local Gibbs relation for a continuous description inside the 
interfacial region, we derive the description of the Gibbs surface in terms of excess densities and fluxes 
along the surface. We introduce overall interfacial resistances and conductances as the coefficients 
' in the force-flux relations for the Gibbs surface. We derive relations between the local resistivities 

^S) \ for the continuous description inside the interfacial region and the overall resistances of the surface 

for transport between the two phases for a mixture. It is shown that interfacial resistances depend 
among other things on the enthalpy profile across the interface. Since this variation is substantial 
the coupling between heat and mass flow across the surface are also substantial. In particular, 
the surface puts up much more resistance to the heat and mass transfer then the homogeneous 
phases over a distance comparable to the thickness of the surface. This is the case not only for 
the pure heat conduction and difi'usion but also for the cross effects like thermal diffusion. For the 
excess fluxes along the surface and the corresponding thermodynamic forces we derive expressions 
for excess conductances as integrals over the local conductivities along the surface. We also show 
' that the curvature of the surface affects only the overall resistances for transport across the surface 

, and not the excess conductivities along the surface. 
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I. INTRODUCTION 



^ . . . 

(~| ' During evaporation and condensation heat and mass are transported through the interface. The common description 

O [ of these phenomena uses assumptions which are questionable. For instance, one usually assumes equilibrium conditions 
^ O ^ ■ at the interface [3| ■ Furthermore one neglects coupling effects between heat and mass transfer . The first one is a 
" zeroth order approximation in deviations from equilibrium and it is desirable to extend it to a more accurate theory. 
y—^ : This is, for instance, done in the monograph by Kjelstrup and Bedeaux (45j . Neglecting coupling coefficients was 
^ ■ shown to violate the second law of thermodynamics for a one-component system (46j . Coupling is important since the 
, corresponding resistances for transport across the surface depend on the enthalpy of vaporization. The significance 
\l ' of this quantity in this context is due to the large difference between the liquid and the vapor values of the enthalpy. 
, Neglecting the coupling coefficients has been shown to imply that the enthalpy of vaporization is zero [4^], which is 
. I clearly incorrect. 

^5 ' The situation becomes even more complicated when one considers evaporation in mixtures. These processes happen 
in such industrial applications as distillation and therefore the precise description is important. Depending on condi- 
tions one can get mass fluxes of components in the same direction or in the direction opposite to the corresponding 
thermodynamic forces. 

There has been done a number of studies of the interfacial transport for one- and two-component systems: ex- 
periments {47-521, molecular dynamic simulations (53l - [57| , kinetic theory [58l - l62j and square gradient continuous 
^ description [63n68|]- AH these works use different approaches, which allows one to investigate different aspects. Most 
5^ , studies were for one-component systems. It was verified that non-equilibrium thermodynamics could be used to 
describe the data by Fang and Ward |47|. lisi , [69| , see also [i^ [s^l- The description of transport through surfaces, 
using non-equilibrium thermodynamics [45| , gives a unifying treatment also for the results of the molecular dynamics 
simulations, kinetic theory and the square gradient theory, both for sirigle component systems and mixtures. An 
alternative approach to interpret the data is the statistical rate theory [t^ [7lj . See for a critical discussion of the 
description of evaporation and condensation the paper by Bond and Struchtrup (72j . 

The results of the molecular dynamics analysis as well as from the square gradient theory give continuous stationary 
state profiles of, for instance, the temperature and the densities tlirough the surface. Both descriptions jssl [sgI [stI, |64. 
Q conclude that the surface is in local equilibrium. This shows that the surface as described by the excess densities 
introduced by Gibbs is a separate thermodynamic system. This makes the systematic description of the surface 
using non-equilibrium thermodynamics given in [45^ possible. One of the points of interest is the dependence of the 
resulting overall interfacial resistances on the continuous profiles. Once we have a description which relates the total 
resistances of the Gibbs surface to the continuous profiles of, in particular, the local resistivities and the enthalpy 
profiles, we can study which aspects of the continuous description are most important.. 
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It is the aim of this paper to obtain the general discrete approach to describe transport properties through and 
along the surface on the basis of the continuous description. As the interfacial region is thin, it is natural to consider 
the whole region as a single entity. In equilibrium, the common method to do this is to introduce the properties of 
the Gibbs surface as excess densities. Away from equililibrium the surface should be described in terms of excess 
densities and fluxes. To obtain these one needs the actual profiles of the thermodynamic densities and fluxes as well 
as the profiles extrapolated from the homogeneous regions. Within such a description, an equilibrium surface may be 
considered as an autonomous phase ^73] ■ Given the validity of local equilibrium it is possible to speak about the Gibbs 
surface also in non-equilibrium states, as was mentioned above. In this paper we will establish the link between the 
non-equilibrium continuous description and the transport properties of the whole surface both through and along the 
surface directly. Crucial in the analysis is that the surface thickness is very small compared to length scales over which 
variables change in the homogeneous phases and along the surface. For fluids this leads to 2-dimensional isotropy of 
the surface discussed in (67| . 

We study the transport through the surface in curvilinear coordinates. This allows us to consider both transport 
through planar interfaces and a phenomenon like nucleation. The curvature of the surface plays an important role in 
nucleation. The agreement of the classical theory [75l - l77j with experiments is in some cases rather poor [t^JsoI] . The 
work in this paper will make it possible to analyze effects due to curvature on the overall resistances for transport of 
heat and mass through the surface of the nucleus in more detail. 

One of the issues this paper focuses on is the overall resistances of the surface to heat and mass transfer through 
the surface. The thermodynamic forces, which lead to heat and mass transfer through the surface, are given by the 
differences between the values of the inverse temperature and chemical potentials divided by the temperature in gas 
and liquid phase. It is known that there are jumps of these quantities across the surface |47h52| . These jumps are given 
in terms of the total heat and mass fluxes using these overall resistances. We refer to a monograph by Kjelstrup et al. 
(45I for a systematic derivation of these linear relations in the context of non-equilibrium thermodynamics of surfaces. 
The other issue studied in this paper is the excess conductances along the surface. Also for these conductances integral 
relations are derived, expressing them in the behavior of the continuous conductivities in the interfacial region. 

The analysis in this paper will only consider stationary states. The resulting overall resistances and conductances 
can also be used to analyze non-stationary behavior. We will neglect viscous contributions to the pressure tensor. 

From the calculation of the overall surface resistances we are able to understand the aspects in the continuous 
description which are most important for the size of the overall resistances. These are found to be the nature of 
the local resistivity profiles and the enthalpy profiles across the surface. Both can be obtained from the particular 
continuous description used. From molecular dynamics simulations there is evidence that local resistivities have a 
peak [5^ [m in the interfacial region. Within the context of the square gradient model studied in (63h68i] this peak 
was modeled by a square gradient contribution. Such a peak clearly contributes significantly to excess resistances. 
The enthalpy is important because the value of the enthalpy is significantly different in gas and liquid phase. The 
overall resistances are found to contain integrals over products of the local resistivity profiles and the enthalpy profiles. 
The location and the nature of the changeover from a large positive enthalpy in the vapor to a large negative enthalpy 
in the liquid relative to the change over and the peak in the local resistivities is important for the size and the sign 
of the overall resistances. The overall conductances are found to contain the excesses of the conductivities for the 
continuous description. 

In Sec. [n] we give a short discussion of the continuous description of the interface using non-equilibrium thermo- 
dynamics. We use the Gibbs relation valid for mixtures in the interfacial region found in the context of the square 
gradient theory, discussed in [gtI, l68j , which we believe to be general. This enables us to obtain the rate of entropy 
production everywhere in the interfacial region. In Sec. |IIIj we define the excess densities needed to describe the 
non-equilibrium Gibbs surface. Furthermore, we show how the small thickness of the interfacial region allows us to 
split all phenomena in separate contributions parallel and perpendicular to the surface in Sec. |IV| . In Sec. jV] we 
derive integral relations for the overall interfacial resistances in terms of the local resistivities and the equilibrium 
enthalpy profiles. In Sec. |VI) we derive the integral relations for the overall interfacial conductances. We further 
derive the corresponding relations for the coefficients for the measurable heat fiuxes in Sec. |VII| . This is done in 
two steps, for normal resistances in Subsec. |VII B| and for parallel conductances in Subsec. |VII C| . In Sec. |VIII| we 
discuss the consequences of the derived relations. In particular we discuss the infiuence of the surface geometry and 
the enthalpy profile across the interfacial region. Finally, in Sec. |IX| we give concluding remarks. 
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II. THERMODYNAMICS OF THE INTERFACIAL REGION 
A. Gibbs relation 

Consider a two phase n-component mixture. Let T be the temperature and p be a scalar pressure. Let hi be the 
partial enthalpy of the z-th component, /i^ be the chemical potential and /im = fii — Hn be the chemical potential 
difference of the i-th component with the 71-th component. Furthermore, let ft., u, s, v be the mass specific enthalpy, 
internal energy, entropy and volume respectively, p = 1/v be the overall mass density and = pi/p he the mass 
fraction of the i-ih component. All these fields depend on the position r and the time t. For a planar interface p{r, t) 
is equal to the parallel pressure P||(r, t). 

The interfacial region is known to have a large density variation over a distance of the order of a nanometer. This 
was the reason for van der Waals [H, HI] to assume, that, for instance, the Helmholtz energy density of the surface 
depends not only on the local properties, but also on their gradients. Due to this dependence the surface tension 
is unequal to zero. For a two-phase system it becomes an additional thermodynamic characteristics and has to be 
accounted in the description of the surface. In the continuous description the role of the surface tension is played by 
the tension tensor 7Q,^(r, t). For a planar surface in equilibrium the only nonzero element of this tension tensor is 
the XX element along the diagonal which equals p±{r) — P||(r). The integral of this difference over the x-coordinate 
normal to the surface gives the surface tension of the planar surface. For a curved surface the tension tensor is no 
longer diagonal. For a fluid-fluid interface it remains symmetric, however. For the relation with the surface tension 
in that case we refer to 68\ . 

A systematic non-equilibrium thermodynamic description requires the Gibbs relation expressing the rate of change 
of the entropy density in those of the other thermodynamic quantities. In [gtI Icsj we used the following form of the 
Gibbs relation 

where v is the barycentric velocity, and d/dt = d/dt + v^d/dxa — d/dt + v-\7 is the substantial time derivative. 
Furthermore we use the summation convention over double Greek indices. Compared to the form of the Gibbs relation 
used in the homogeneous regions (83 |. it has an additional contribution proportional to the divergence of the tension 
tensor 7q.^. This contribution is only nonzero in the interfacial region. The use of it emphasizes the difference 
between the surface and the homogeneous regions. The above expression is equivalent to the one given in j63j for the 
one-component two phase fluid. In [g^, HI] we motivate the choice of the new term in the interfacial region. This 
motivation is not a derivation, however, and the validity of the above Gibbs relation is therefore postulated. The 
justification of the specific choice of the Gibbs relation is on the one hand a comparison with experiments and on 
the other hand the thermodynamic consistency of the results. Such a consistancy was found to be lacking for other 
choices. 

The next step in the thermodynamic description is to give the relation between the thermodynamic potentials. In 
a bulk region they are known to be homogeneous functions of the first order. It is not obvious that this property 
remains valid inside the interfacial region. It follows however from Eq. (jll.l[) that the internal energy density has to 
be a homogeneous function of the first order of the entropy density, volume and mass densities. Indeed, all the terms 
in Eq. pi.ip but the last one have the common form. The last term is proportional to the gradient of the tension 
tensor, which is not an extensive quantity, however, so it does not contribute to the internal energy. We therefore 
may write the following relation both for the homogeneous as well for the interfacial regions 

n-1 

u{r, t) = ^„(r, t) + ^ M^n(r, t) Ut, t) - p{r, t) v{r, t) + r(r, t) s{y, t) (II.2) 

i=l 

We emphasize the role of the interfacial region here. While the thermodynamic potentials are related via the stan- 
dard thermodynamic relations, the Gibbs relation in the interfacial region is different from the bulk Gibbs relation. 
Important is that the standard thermodynamic relations contain p[r, t), which is equal to the parallel pressure for a 
planar interface and which is very different in the interfacial region from the asymptotic values away from the surface. 
Substituting Eq. (|II.2p into Eq. (III.ll) we obtain 



(11.3) 
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This is the Gibbs-Duhem equation for a two-phase multi-component mixture. Note that since fj,in = fJ^i — fJ.n and 
J27=i 6 = 1 we have 

, n—l J ^ J 

dt dt dt ^^^-^^ 

which is the usual contribution to the Gibbs-Duhem equation associated with the chemical potentials. To simplify 
the notation we will not usually indicate the r, t dependence. 

For a stationary state the derivative d/dt gives zero and Eq. pi.3|) takes the following form 

y dxp dxp ^ * dxp " dxa j 

where aap = pSap + "jap is the thermodynamic pressure tensor. 

The equation of motion for a stationary state takes the following form: 

dvp duap a(g-r) 
pv„— ^ = p—^ (11.6) 

OXa OXa OXj3 

where g-r is the gravitational potential. Contracting this equation with v^ and dividing by p the result may be 
written as: 

Using it together with Eq. (|II.5P we get 



v-Ee.fv^-/i,vl)=0 (IL 



which is the Gibbs-Duhem equation valid both in the homogeneous regions as well as in the interfacial region in a 
stationary state. Here pi = pi + v^/2 — g- r and hi = hi + — g-r. Note, that hi = 'jli + Tsi, where Si is the partial 
entropy. 

B. Entropy production 

The entropy balance equation is 

P^--V-J,+a, (II.9) 
dt 

with the entropy flux = Js.tot ~ psv and the entropy production as- Given the Gibbs relation Eq. pi.l[) it is shown 
in (stI . l68j that the entropy flux and the entropy production are given by the following expressions 

Js = ^ ^Jg - ^ Pin Jj^ (II. 10a) 

n-l 



where 3q and are the heat and diffusion fluxes 



Jq = Jg - pve -pv = Je - Jm Qi + v^/2 - g-r) 

(11.11) 

Ji = Pi (Vi - v) = J^, - ^i Jm 
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The diffusive mass fluxes satisfy the relation Y^^=i — 0- In turn, the energy flux Jg and the mass fluxes J^. = pi 
and = p V are convenient quantities since in stationary states 

V-Je = 0, V-J^, =0, V-J,„ = (11.12) 

Furthermore, it follows from Eq. (III.9I) that in a stationary state 

Gs = V-is,tot = V-J, +pv-Vs = V-(J, +psv) (11.13) 

Using the Gibbs-Duhem equation (jll.5|) and the conservation laws, under stationary state conditions, it is possible to 
show that 



1 " ~ 



The fluid-fluid interface has a two-dimensional isotropy. A detailed discussion of what this means in a continuous 
3-dimensional description is given in [g^] . The essential element in this discussion is the fact that the interfacial region 
is thin compared to the radius of curvature, 5 <^ ^. This implies that the 3-dimensional vectorial fluxes and forces 
in the entropy production (jII.14[) split up into 2-dimensional vectorial components parallel to the surface and scalar 
components normal to the surface. We may therefore write 

as = CTsji + (Ts,i (11.15) 

where the parallel and normal contributions are 

^ n ~ 

CTs.ll = Je,|| • V|| - - ^ J?„ll • V|| Y (11-16) 

2 = 1 

n ~ 
i=l 

Two-dimensional isotropy now implies that the parallel and the normal forces and fluxes do not couple and the force- 
flux relations for parallel and perpendicular contribution are independent. As we will see below, it is convenient to 
write these relations in terms of conductivities for parallel contributions and in terms of resistivities for perpendicular 
contributions. We therefore have the following force-flux relations in the parallel direction 



1 " f Jl- 



(11.18) 



and in the normal direction 



1 " 

(11.19) 



4=1 



The 2-dimensional isotropy implies that transport coefficients for both contributions (conductivities for parallel contri- 
bution as well as resistivities for normal contribution) are scalar. The off-diagonal coefficients of both sets of force-flux 
equations satisfy the Onsager reciprocal relations. Inside the interfacial region the normal resistivities are different 
from the corresponding parallel resistivities. Away from the interfacial region the fluid has a 3-dimensional isotropy 
and the corresponding resistivities (or conductivities) for the parallel and the normal direction become equal. In 
order to be consistent with the sign conventions used in the monograph by Kjelstrup et al. [45| we changed the sign 
convention we used in (gtI . Issj . For this purpose we changed the sign of the chemical forces which results in a change 
of the sign of the cross coefficients between the energy flux and the mass fluxes. 

We emphasize the reader's attention on the following difference in the convention. It is common to give the 
constitutive force-flux relations for n-component fluid in terms of n independent fluxes, e.g. the heat flux Jg and 
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n — 1 diffusive fluxes . Those fluxes are linearly independent and the corresponding matrix of transport coefficients 
(resistivities or conductivities) is non-singular. If one uses n-\-l fluxes or forces as variables, the corresponding matrix 
of the transport coefficients has n + 1 rows and columns and is therefore singular. This makes it impossible to convert 
resistivities into conductivities that easily, so the corresponding representation should be fixed. The use of particular 
representation (resistivities or conductivities) is dictated by the surface symmetry and will be apparent later. The 
singularity of the matrix of transport coefficients does not however restrict us much. For instance, there is one-to- 
one correspondence between the coefficients of n-squared matrix and n + 1-squared matrix (see e.g. Appendix [C] ). 
Furthermore, the singularity of the transport coefficients matrix does not affect the Onsager reciprocal relations, cf. 
[S^ . ch.VI, §3]. We can therefore operate these coefficients normally, unless we want to invert the matrix. 



III. THE GIBBS SURFACE 

The Gibbs surface is described with the help of excesses of different quantities in the interfacial region. In equilibrium 
only excesses of the densities are used in the description. Non-equilibrium involves many other quantities which vary 
through the interfacial region and which therefore require a proper definition of the corresponding excess quantity. 
In this section we will define the excess densities in non-equilibrium. In the following sections the definition of the 
excess resistance, the excess conductance, and the excess fluxes will be given. 

Because of the variety in nature of the quantities we want take the excess of, we shall use the unified notation for 
the excess. Namely, we will use ^riq] to denote the excess of a profile q{r), where the subscript r indicates the nature 
of the the quantity q. For example, £y will indicate the excess of a density and £j will indicate the excess of a flux. 
Even though both <Sv and €j represent the property of the Gibbs surface, they are defined differently. 

The definition of an excess requires the normal direction n to be defined in the interfacial region. The surface may be 
curved and we may introduce curvilinear orthogonal coordinates r = {xi,X2,X3) with xi being the normal coordinate 
and r|| = (x2, X3) being the tangential coordinates. For the stationary states these coordinates are independent of the 
time. 

Let be a variable defined in the interfacial and the homogeneous regions. Furthermore, for any variable (j>, let (j>^, 
where superscript b stands either for £ or for g, be the variable 4> extrapolated from the bulk to the surface region. 
The extrapolation is done using the description in homogeneous phases. Outside of the interfacial region (j)^ and (p are 
almost identical but inside the surface (j>^ in general differs from (f>. In order to make this comparison more precise, 
let a;^''*(r||) and a;^''*(r||) be the boundaries of the interfacial region at the gas and liquid side respectively. These 
boundaries are chosen such that the extrapolated value cffi differs some small fraction (like for instance 0.1 promille 
or 0.01 promille) from the actual value (j) along these boundaries. Given the small size of this difference we can use 
the following identity for the extrapolated variables 

0''(x'''«(r„),r,|) = 0(x''^^(r„),r|,) (III.l) 

The location of the boundaries of the interfacial region depend on the variable and the accuracy used. For 
a sensible choice of the variable and the accuracy, all choices are equivalent. In contrast to equilibrium, in non- 
equilibrium description all the bulk variables vary in space and so do the extrapolated quantities. This is true both 
for densities, like the density p, which vary only slightly in the bulk phases, and for intensive quantities, like the 
temperature, which may vary significantly even in the bulk phases. If one extrapolates any function of a number 
of variables the result equals the original function of the extrapolated variables: 

F''(...,0,...) = F(.. .,/,...) (III.2) 

This is a consequence of the fact that the equality is obvious away from the surface. For instance, (/ij/T)^ = /if /T^. 

We note however, that even though we introduce the extrapolated variables in such a way that Eq. pil.ip and 
Eq. (jIII.2[) are true, the use of any numerical procedure to obtain the extrapolated variables will introduce small 
errors. This is true in particular for non-equilibrium when the extrapolated variables are not constant. This is 
however a technical question, not a fundamental one. Usually the extrapolation involves polynomials in order to fit 
an actual curve, which introduces a non-zero error in the extrapolated curve. 

On the basis of some criterium one chooses a dividing surface in the interfacial region. The location along the xi 
coordinate is denoted by a;''(r||). This location is usually, though not necessarily, between x^'^{r^^) and a;^'*(r||). We 
will discuss below which criteria are used for this choice. 



7 



A. Excess density 

If ^(r) is a density per unit of volume in the S-dimensional space, we can define its excess €y [</>] as (8 Si] 

£y[0](x^(r||),r||)^-l- / dxif)i[)2[)3<A^"(r;x^(r,|))) (IIL3) 

where 

0-(r;x^(r,|)) ^ <t>{v) - ^^(r) e(^^(r,|) - x,) - 0^(r) e(xi - a:^(r||)) (III.4) 

Furthermore, \]a = ^aixi^Y^^) are Lame coefficients for the curvihnear coordinates and = \]a{x'^ {y\\)tY\\) for 

1 = 2,3 are Lame coefficients for the curvihnear coordinates along the dividing surface. Given that 4>{v) is a density 
per unit of volume, the excess [4'] i^'' is a- density per unit of surface. The excess depends on the position 
of the dividing surface a;''(r||). For ease of notation we will not further indicate the dependence of a;*, a;^'" and x^'" on 
Ty explicitly. 

The key quantity in the non-equilibrium description is the entropy production. As for any density, one can introduce 
the excess entropy production for the Gibbs surface. It will determine the dissipation rate of the whole surface. 
Consider the entropy production given in Eq. pi. 141) . All the terms have the form J-Vro, where J is a flux and m is 
some scalar function. According to Eq. (|II.12p V-J = in the stationary states, so J-Vro — V-{3w). We show in 
Appendix [X] that this leads to 

£y [J-Vnj] = {J^vjY - {J^vjf + £y [V||-(J||n7)] (IIL5) 

where all the functions on the right hand side are evaluated at r*^. We see, that the excess entropy production of the 
Gibbs surface in the stationary states can be split into two contributions: the one takes into account only perturbations 
perpendicular to the surface, the other depends only on the perturbations parallel to the surface: 

£y [a,] = (£y ] + £y [(7,,^ ] (IIL6) 

We note, that such a separation does not require the surface to be two-dimensionally isotropic. 



IV. A THIN INTERFACE 



Substantial simplification in the expression for the excess entro py production can be made if one consider typical 
conditions, when the interfacial region is thin. As was discussed in |67| . the interfacial region breaks the 3-dimensional 
isotropy of the system. In addition to a typical macroscopic size of the problem D, there exists the microscopic size 
d <^ D, the surface width, which is of the order of few nanometers. 

First, we shall make a note on the definition of the excess density. For a thin surface the Lame coefficients 
f)a = f)a(a;±,r||) do not vary much along x^_ coordinate. To a relative order {d/D) they are equal to 

f)„(x,,r||) = [,„(x^r||) + d£^ + ... (IV.l) 

Substituting this expansion into Eq. (|III.3[) . we get 

€y[<^]-£[0]+O(d/i?) (IV.2) 

where 

(£[0]=/ dxif)i<?!)'=^(r;x') (IV.3) 



For each flux in Ec^. pi.l2p we can write Vj^ Jj^ + V|| • J|| =0. This gives an approximate relation for the order of 
magnitude 

\A±IA ^ \hIA (IV.4) 
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where and denote a typical change in the corresponding direction. There are quantities which change drastically 
over the distances of the order d in the direction perpendicular to the surface, so A^Xj^ ~ d. However, a significant 
change of any quantity along the surface may happen only on a length scale D, which is of the order of either the 
radii of curvature or the system size, so A^x^ ~ D. Because of this property of a surface, we may not expect that 
[Aj^J^I ^ |A|| J|||, i.e. that the change of the parallel component of a flux on a macroscopic scale along the surface is 
much larger then the change of the perpendicular component of that fiux on a microscopic scale through the surface. 
For the fluxes for which changes | Aj^ | and | Ay Jy | are of the same order of magnitude A J, Eq. (IIV.4P takes the form 
AJ/d w AJ/D, which can hold only if A J = 0, since d <^ D. This means that both A^J^ — and Ay Jy = 0. If 
Aj^J^I |A|| Jyl this statement becomes even stronger. We then may require that for a thin surfac43 



V,rJ„ r 





: 



(IV.5a) 
(IV.5b) 



Thus, for a surface with a thickness much smaller than the radii of curvature, the stationary state condition V-J = 
has the form of Eq. (jIV.5[) in an interfacial region. 
The extrapolated fluxes J** satisfy the same equations 



V 



= 
= 



(IV.6a) 
(IV.6b) 



since the extrapolated flux flelds also satisfy V- = 0. 

Consider now Eq. pV.Sap and Eq. (IIV.6al) . Both of them are first order ordinary differential equations, the solutions 
of which contains additive constants. These constants must be the same, since according to Eq. (jllLTI) J^(r^'") = 
J^{y^'^) at the boundary points. It follows that and are the same functions: 

Jf(r) = Jf(r) = J,(r) (IV.7) 

Note that Eq. pV.7|) does not lead to the relation Jj^(a;^'*, ry) = Jj^(j;^'*, ry). Eq. ljIV.7p is the relation between values 
of different functions at the same point but not the relation between values of the same function at different points. 
However, it follows from Eq. (jIV.Sal) that in curvilinear coordinates d{\)2 f)3 J±)/{dx^) = and therefore 

i}2ix^) U^^) = const = ()^ f,^ J(x^) (IV.8) 

but not = const. 

Consider now Eq. pV.Sbp and Eq. pv.6bp . It follows that 



V||-(J||G7) =J||-V||tn 
Vy(jM 



(IV.9) 



Applying both Eq. pVJl) and Eq. ffPT^ to Eq. (111131) we get 

€ [J-Vw] ^ J^{vj^ - ccjS) + (£ [ J|| -Viinj] (IV.IO) 

Applying Eq. pV.lOp to each term in Eq. (III.14P for the entropy production we obtain the general form of the excess 
entropy production for a surface in stationary state 

1 



El 



T3 



(IV.ll) 



'e, II 



•V||i 



'6, 



• V ^ 



We see, that for a thin surface, the excess entropy production splits into independent contributions which are caused 
by forces perpendicular and parallel to the surface. This property has the same origin as the 2-dimensional isotropy 
of the surface j67|. 

It is desired to have the excess entropy production in the form Q; [cts ] — Jj Xj, i.e. as a sum of independent 
contributions, each of which can is a product of a flux and a force. However, as one can see from Eq. pV.lip the 
parallel contributions do not have such a form. Further information about the symmetry of a system is required. We 
will therefore consider the cases when the forces are applied perpendicular to the surface and parallel to the surface 
separately. 



^ For the special case of a system with planar surface in cartesian coordinated with all the fluxes directed perpendicular to the surface, 
these equations follow straightforwardly. 
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V. INTEGRAL RELATIONS FOR TRANSPORT THROUGH THE SURFACE 

Consider forces applied only in the direction perpendicular to the surface. Then CTsji = and the local entropy 
production cr<, = as,± which is given by Eq. pi.l7|) . The corresponding force-flux relation are given by Eq. (IIL19|) . 
The excess entropy production given by Eq. (jIV.ll[) simplifies to the following 

[ ] = ^ - ^) - E ^e. . (I - §) (v.i) 

The fluxes in Eq. (|V.1|) are evaluated at the dividing surface x'^ . Furthermore, = T^{x^) and = T^{x^) are the 
temperatures extrapolated from the liquid and gas to the dividing surface. The analogous meaning have /If and /I^. 

Following the traditional procedure in non-equilibrium thermodynamics, we write the force-flux relations based on 
the form of the excess entropy production (jV.l[) as: 
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(V.2) 



We now compare Eq. (jV.2p with Eq. pi.l9[) . Eq. (|V.2I) for the force- flux relations has the following form 

w\x')-vd3{x')^ Rl{x')J^{x') (V.3) 

The resistance i?'^ in general depends on the choice of the dividing surface. Eq. (III.19P has the following form 

V^vj{x)^rl{x)J^{x) (V.4) 

Let us apply £ to both sides of Eq. (jV.4|) . We note, that we use the excess operator £, given by Eq. (jIV.3|) exactly, 
not any other operator which is equal to € only to the order d/D. Using Eq. (jIV.8|) and the fact that in curvilinear 
coordinates Vj_ = d/dxi, one obtains 

m\x') - m^ix') = (Br [r^ ] (x") J^x') (V.5) 

where we introduced the excess resistance operator (£r. For a resistivity profile r(r) it is defined as 

£fi[r]^[,^(,^ / dx^-^r^^{r;xn (V.6) 

Jx^^-" 1)21)3 

Comparing Eq. (IV. 51) with Eq. (|V.3I) we therefore conclude that 

Rl^€n[rl] (V.7) 

This is the general form of the resistivity integral relation for the transport normal to the surface. 

Given that the forces are applied only in the direction perpendicular to the surface, Eq. (|V.7p is valid for the 
interfacial region of any thickness. We now make use of the fact that the surface is thin. For such a case we expand 
the Lame coefficients l)a{x) around and obtain 

iBR[r]^€[r] + Oid/D) (V.8) 

so that to a relative order {d/D) the excess resistance [r ] is given by £[''']. 

Applying the above procedure to each term in Eq. (jlLlQp and Eq. (|V.2p we may conclude that 

^'-qi.-L ^ L qt,l- J (V.9) 



to the order {d/D). Eq. (jV.Sp gives integral relations for the resistivity coefficients for the case that one uses the 
total energy flux to describe the transport through the surface. They show that resistivities normal to the surface are 
additive, as one would expect. 



10 



VI. INTEGRAL RELATIONS FOR TRANSPORT ALONG THE SURFACE 

Consider now the perturbations applied only in the direction parallel to the surface. Then as^± ~ and the 
local entropy production ag — (jg y which is given by Eq. pL16[) . The corresponding force-flux relation are given by 
Eq. (Iimi) . 

A. Excess entropy production 

The excess entropy production given by Eq. (jIV.lip simplifies to the following 



Je, II • V|| — 



T Y7 



(VI.l) 







1 


/ 










1 


/ 







This form of the excess entropy production does not allow us to write the constitutive force flux relations yet. 
However, the fact that there are no perturbations in the direction perpendicular to the surface, allows us to simplify 
this expression further. 

Each term in Eq. (IVI.ip has a form € [ Jy ■ Vytu], where the possible perturbation is controlled by a scalar function 
Tu{xj^, r||). For perturbations applied only in the direction parallel to the surface, ■uj is independent of x^. The same 
argument is applicable to the extrapolated property w''. We show in Appendix [B] that this leads to the following 
relation 

€y[J|,-V||n7] = (£j[J||]-V||n7 (VI.2) 

where we have introduced the excess flux operator €,/. For a parallel component Jy = (J2, J3) of a flux J = (Jj^, Jy), 
the excess parallel flux is defined as (£,/ [ J,, ] = {^j^ [J2] , Sja [J3]), where 



(VI.3) 



where 

J^^(r; x"") = J(r) - J^r) e{x' - xi) - J^(r) 6(2:1 - x') (VI.4) 

As discussed in more detail in possible excess fluxes normal to the surface, calculated in a reference frame that 
moves along with the surface in the normal direction, play no role in the discreet description of transport through the 
surface. Here we note that for the special case of stationary states the surface is not moving in the normal direction. 
For a thin surface one can expand the Lame coefficients around x'^ and obtain to a relative order (d/D) that 

€j[J,] = €[J,]+Oid/D) (VI.5) 

where operator £ for a flux Jy is still given by Eq. pV.3l) . 

Following [sF] we denote the excess parallel flux €j [ Ji, 1 as JjJ. 

Consider now V||ti7 which appears in Eq. (jVI.2p . In [68|] we discussed how to deflne the temperature T** and the 
chemical potentials /if of the surface in terms of the surface tension and the relative adsorptions. These quantities 
are independent of the choice of the dividing surface when the surface is planar. For a curved surface differences of 
the relative order of the distance between the dividing surfaces divided by the curvature {d/D) may appear. The 
quantity uj considered above can be written as a function of the temperature and the chemical potentials alone. We 
can therefore evaluate this quantity at the temperature T" and the chemical potentials /if of the surface. This then 
gives ccj^, which therefore will be independent of the choice of the dividing surface to a relative order d/D. As we 
have discussed, the variation of vj along the surface is the same as the variation of vj^ along the surface. It therefore 
follows that the variation of tn* along the surface should also be the same to this order. Eq. (jVI.2p takes the following 
form 

£[J„-V„n7]=J^V„t^^ (VI.6) 

to this order. 
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Furthermore, we identify the excess entropy production in Eq. (|VI.ip as the surface entropy production: (t| h . With 
this notation, substituting Eq. (|VI.6p into Eq. (|VL1I) and taking into account the terms to the order d/D we get for 
the excess entropy production 

<„=J:,,-v„--5:j|,,,,v„^ (VI.7) 

which has the form used in [86i] . 



B. Integral relations 



According to the traditional procedure in non-equihbrium thermodynamics, the force-flux relations which follow 
from the excess entropy production (jVL7|) are 



i=l 



1 " 

+ ^ ^ij, II 



V|| j,^ 



(VI.8) 



Due to the isotropy of the surface the conductance coefhcients LJ| are scalars. 

We now compare Eq. (|VI.8I) with Eq. (III.18|) . Eq. (IVI.8P for the force- flux relations have the following form 

3;ixn = mx^)y,^' (VI.9) 

Both Ljj and £ [ Jy ] in general depend on the choice of the dividing surface, while VyHj* does not. Eq. pi.l8l) has the 
following form 



J|i(x) =£l{x)V^^m{x) 



(VI.IO) 



Let us apply €j to the both sides of Eq. (|VI.10I) . We note, that even though the excess operator (£j given by 
Eq. (jVI.3|) is equal to the operator (B to the order d/D, we will first keep the higher order terms. Direct substitution 
yields 

3lixn = eL[il]{xn-^,w^ (VI.ll) 

where we introduced the excess conductance operator (Bl- For a conductivity profile £{r) it is defined as a matrix 

n \ 



/M 
V 







t)2 



77 / rfa;if)i7^ 



(VI.12) 



f)3 



Next we would like to compare Eq. (IVI.9P with Eq. (|VI.11[) in order to relate the interfacial conductance L|[ to the 
excess conductance (Bl [ ] ■ This however is not possible directly, since in general the former one is a scalar while 
the latter one is a matrix. However, if we take into account that the surface is thin, we expand the Lame coefficients 
i)a{x) around x^ and obtain 



_ f€[£] 



'^^+Oid/D) 



(VI. 13) 



so that €l [£] = to a relative order {d/D). We note, that to this order the excess [i] becomes scalar. 

Comparing then Eq. (jVI.ll[) with Eq. (|VI.9p we can conclude that 



(VI.14) 



to the order {d/D). Neglecting curvature effects, this is the general form of the conductivity integral relation for the 
transport parallel to the surface. 
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Comparing Eq. pTTS)) with Eq. dVLSl) and using Eq. (jVI.14p we may conclude that 





= (£ 




T ^ 

qi,\\ 




\ 

I q-i;\ 


T <^ 


= £ 





(VI.15) 



Eq. (jV.8[) represents integral relations for the conductivity coefficients for the case that one uses the total energy flux 
to describe the transport along the surface. They show that conductivities parallel to the surface are additive, as one 
would expect. 



VII. MEASURABLE HEAT FLUXES 
A. Local transport coefficients 

It is convenient to write the local entropy production also in terms of the measurable heat flux , which is defined 

as 

n n 

j;=Jg-^M«=Je-^^^J«. (VII.l) 

i=l i=l 

In terms of measurable heat fluxes the expression for the entropy production becomes 

^.=J^V--^J,. V^-/i„V- (VII.2) 

where /ij„ = hi — /i„. Eq. pi.lSp . Eq. (jll.141) and Eq. (jVII.2p . which give the entropy production in stationary 
states, are equivalent expressions. To second order in the deviation from equilibrium one may, and we will, use the 
equilibrium values of hin in Eq. ()VII.2p . The resulting expression for the entropy production is then only equivalent 
to the other two to this order. 

As in Sec. |II B| we can split the local entropy production into parallel and perpendicular contributions based on 
the two-dimensional isotropy of the interfacial region, see Eq. pi. 151) , so that 



(^SA = Jq,||-V||^ - ^^^,r\^\\^ - ^m.egVii^ j (VII.3) 

= '^.UV4-E'^^- fv,^-/.„,e,v4) (VII.4) 



1=1 

n-1 

^ S.J- " (T I • _L rji / J " £, _L I • _L rji ' "t /t , et/ • _l rjn 

i=l ^ 

This leads to the following force-flux relations in the parallel direction 

n-l 



i=i ^ ' 



n-l 



(VII.5) 



and in the normal direction 



rp ' ^ ^ IJ \ • W rp ' --tti^t^i^ II rp 

1=1 



^ n — 1 



i=l 
n-l 



(VII.6) 
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Further on, for the transport through the surface we wih use only resistivities and for the transport along the surface 
we will use only conductivities. We have therefore omitted subscripts ± and ii for corresponding transport coefficients 
to simplify the notation. 

We now relate the local resistivities associated with the measurable heat flux to the local resistivities associated 
with the total energy flux. Comparing Eq. (jVII.Gp with Eq. (jII.19[) we obtain 



^qi ^qq ^i.eq^ 



i = l,n 



(VII.7) 



"^ji — T"qi hj Qq Tjq hi ^q -\- Tqq llj ^q lli ^q^ i^j — 1, ?T- 



where rgo, Tg^, and defined as 



' 99 — ' qq 



n-1 



Tqi = < 



fc=i 

n-l 



fc=l 
/ n — 1 71 — 1 



n-1 



k=l 1=1 k=l 
n—l n—1 n— 1 

fe=l 1=1 k=l 
n— — 1 n — l 

Y^rl^i ii-^^k rli, 



(VII.8 



i = n, j = l,n — l 



k=l 1=1 
n—l n—l 



k=l 



i = l,n — l, j = n 



I, J = n 



k=l 1=1 



We now relate the local conductivities. Comparing Eq. pi.l8|) and Eq. (jVILSP we get 



n-l 



99 ~ 99 



n-l 



j=l i=l 



^qi ^ ^gi ^ ^ ] ^jn,eq ^ij i 



i — l,n 
i,j = l,n 



(VII.9) 



^ The details of this procedure are given in Appendix [U] . 
^ The details of this procedure are given in Appendix [D] 
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where Iqq, £qi, and £ji are defined as 



f = f 



^qi, i = l,n-l 



k=l 



(VII. 10) 



i = n, j — 1, n~l 

-E j = j = n 

fc=i 

n— 1 n— 1 

E E ^fc" = " 



fc=l 
n-1 



k=l 1=1 



In the rest of the section we will derive integral relations for the measurable resistivities and conductivities. As 
such, we will relate the measurable interfacial resistances to the measurable resistivity profiles and the measurable 
interfacial conductance to the measurable conductivity profiles. The reason for this is that the coefficients which 
correspond to the measurable heat flux can be obtained experimentally. 



B. Transport through the surface 

The excess entropy production due to transport through the surface can be written in terms of measurable heat 
fluxes extrapolated from either gas of the liquid phase tho the surface. The extrapolated measurable heat flux is 
defined similarly to Eq. (IVII.1[) : 



•J q^l_ — "Je 



(VII.ll) 



In equilibrium the fluxes are equal to zero. To linear order in these fluxes we may replace the enthalpies by their 
equilibrium values hi^ eq = hi^eq — g'l"^. In the following equations one should therefore use hi^ eq in stead of hi. This is 
similar to the use of the equilibrium enthalpy difference profiles in Eq. (|VII.2p . Using Eq. (|VII.11[) . the excess entropy 
production (|V.ip becomes to linear order 



J', a 



1 



1 

T9 



i=l 



T9 



"•i.eq 



1 



1 

T9 



(VII.12) 



The measurable heat flux Jq]l is the one calculated for the gas side of the surface. The similar expression can be 



written for the measurable flux V calculated for the liquid side of the surface. 
The corresponding force-flux equations for the gas-side are 



1 L ^ K"9 j'-9 

rpi rpg ^^qq "^9.-L 













T9 1 




T9 J 



Comparing Eq. (|V?2|) with Eq. (jVILlSp we find 



E 

n 

= ^;^^:! + E 

i=l 



Rqi Jii 



(VII.13) 



?9 



"l, eq ^'-qq 



(VII. 14) 



3, eq "qi 



I, eg ],eq "qq 
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Using integral relations Eq. (jV.9|) for absolute resistivities together with Eq. (jVII.14[) and Eq. (IVIL7p . we find 
integral relations for the measurable resistivities: 



^qq ~ ^ [ ^qq ] 

Kf =^[rq^] + <^\ihlg-h^. 



(VII.15) 



(M.eq K,eq) 1~jq 



He, 



Keq-hj-ei)iKeq 



For the record we also give integral relations for the measurable interfacial resistances on the liquid side of the surface: 

Rqq — 2 [ fqq ] 



i,eq hi^eq) Tqq 



(VII. 16) 



=(S[r,,] + (S ihl,^-h,,eq)r,q 



(^j,eq ^j.eq) fqi 



gg 



Note, that the excess resistances are linear in cc*. For example, if we consider the two different positions of the 
dividing surface xf and x^, then for the excess resistance to the heat transfer we get 



RqqK) ~ Rqqi^2) ^ K - A){rl^ - T^^) 

The similar relations hold for other resistances. 



(VII.17) 



C. Transport along the interface 

In contrast to the transport perpendicular to the surface, in the case of the transport parallel to the surface we are 
not interested in the extrapolated quantities. The reason for this is, as was mentioned above, that the surface is a 
separate entity and the transport along the surface cannot be determined by the properties of adjacent phases. This 
observation is in agreement with the hypothesis of local equilibrium of the surface, which we verified in [68[. In view 
of this the measurable heat flux along the surface should be defined differently from the one normal to the surface. 
The measurable heat flux is the total energy flux corrected for the comoving energy due to the mass flux. As the 
surface is a separate entity, we should substract the excess enthalpy flux of the whole surface from the excess of the 
total energy flux. 

The enthalpy of the surface is defined using the excess densities. If 0(r) is a specific quantity per unit of mass, the 
corresponding surface specific quantity is defined as [86| 



where p is the mass density. Given that 0(r) has a dimensionality of a quantity per unit of mass, the surface quantity 
has a dimensionality of a quantity per unit of mass as well. 
We define the measurable heat flux along the interface as 

n 

J^'.lf ^Je,„-E^^'J|..ll (VII.19) 
i=l 

where the surface enthalpies are defined using Eq. (jVILlSp . We emphasize the difference between the perpendicular 
and the parallel transport here. In the case of the transport normal to the surface, we have the measurable fluxes on 
the gas side of the surface, Jq^l, and on the liquid side of the surface, Jq]i^ which are different from each other. In 
the case of the transport along the surface, we only have the measurable flux J^'if. 

In terms of the measurable heat flux the excess entropy production due to transport along the surface (|VI.7[) can 
be written using Eq. (|VII.19p to linear order: 



V h" V — 



(VII.20) 
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where we replaced the surface enthalpies hf with their equilibrium values hf which is correct to linear order. 



The corresponding force-flux relations are 

1 " 



qi 



\7 — -\-h^ V — 



(VII.21) 



where we expressed fluxes in terms of forces, as we did it before for parallel contributions. Comparing Eq. (jVI.8[) and 
Eq. (jVII.21[) we get for the measurable conductances 



T Is 

^qq 


- ^qq 










^qi 


— ^qi 


n 















1=1 J=l 



(VII.22) 



J2 



Using Eq. (IVI.15|) for absolute conductivities together with Eq. (|VII.22|) and Eq. (|VII.9|) . we find integral relations 
for measurable conductivities 



71 

L'qq ^^[lqq]~Y.^[(hl,,-\eq)^ 
j = l 

n 

L'q! =^\^1^^]-T.'^[Keq-H.eq)i 



vjq -T- ^qj 



n n 

) ] + E E ^ [ - heq) (Keq - heq) 



■qq 



(VII.23) 



Ljf depends linearly on the position of a;'', i^f depends on a;'' quadratically, and L^^ depends on cubically. 



VIII. DISCUSSION 



A. Dependence on the enthalpy profile 

An important property to analyze is the dependence of the thermodynamic quantities on the reference state. All 
the thermodynamic potentials are defined with respect to some reference state. The behavior of the resistances 
and resistivities r"^ on one hand and R'^ and r' on the other hand is different in this respect. Both R"^ and r*^ 
coefficients depend on the reference state, as they are associated with the absolute fluxes. The coefficients R'^ and 
r' are associated with the measurable fluxes and therefore independent of the reference state. The similar arguments 
are valid for the interfacial conductances and conductivities. 

Consider integral relations for the measurable resistances, Eq. (jVII.lSp or Eq. (|VII.16p . Each component of the 
resistance matrix R'^ contains the excess of the local resistivity profiles coupled to the enthalpy profiles. Since all 
the terms on the right hand side of these equations are evaluated at equilibrium, the resistances depend linearly on 
the position of the dividing surface. 

The first term in the expression for every resistance is the excess of only the corresponding resistivity profile. As we 
discussed in [ij, each resistivity profile contains a peak due to the interfacial contribution [sij. It is the size of this 
peak which controls the magnitude of the first term and makes it positive. This contribution is the only contribution 
to the interfacial resistance R^'^^ = R'^'^ to heat fiow. In this sense, the resistance of the surface to heat flow is similar 
to the heat resistance of a finite bulk layer. 

The resistances to the mass transfer due to the temperature difference, contain the excess of an additional 

contributions, however. This contribution depend on the variation of the enthalpy profile across the interface /if — 
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(a) 



(b) 



FIG. 1. Schematic representation of the enthalpy profile hi{x) (solid line) and the resistivity profile rqq(x) (dash-dotted line): 



1(a) the resistivity profile is shifted to the left; 1(b) the resistivity profile is shifted to the right 



hi eq and the local heat resistivity profile. The variation of the enthalpy difference changes from zero on the gas side 
of the surface to the value of the enthalpy of vaporization Aft,""^ = h\ — on the liquid side. This gives an enormous 
contribution, negative for R'^^^ and positive for R^j^ , to the resistances to the mass transfer due to the temperature 
difference. The magnitude of this contribution depends on the relative position of the enthalpy and the resistivity 
profiles. This is illustrated in Fig. [T], where the peak in Tqq hardly contributes to the integral of the product in 
Fig. [l(a)| , while it contributes substantially in Fig. [l(b)| . Depending on the sign of the coupling coefficient, the 
other contributions may be positive or negative. 

In contrast to the resistances, the conductances are not linear with respect to the position of the dividing surface. 
This is due to the dependence on the excess enthalpy of the surface on the position of the dividing surface. The other 
arguments given above for the resistances, are similarly applicable to the conductances. 



B. Influence of the system geometry 

Here we consider some typical cases of the surface geometry. As the analysis above is general, we apply it to planar 
interfaces and droplets. 

For planar interface the lame coefficients fli = 1)2 = f)3 = 1 and the excesses defined above in Eq. (jIII.3p . Eq. (|V.6p . 
and Eq. (IVI.12I) take the following form: 

^R[r]{x') ^^{r}{x') (VIII.l) 

where 

'P{g}= / dxq''''{x]x') (VIII.2) 

is the excess operator for planar interface and q^'^ is defined similarly to Eq. pil.4[) . We see, that for planar surface 
the approximation made throughout the paper is satisfied exactly. 

For droplets, which have spherical symmetry, we have f)i = 1, [}2 = x, [}3 = xsmO where x is the position coordinate 
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and 9 is the polar coordinate. Then, the excesses defined above take the following form: 



rX' /^s\2 

(Bii[r]{x') ^ j "^^(^1 ^'"(^;^') (VIII.3) 
= / dxr'=(x;x') 



J a;" 

where x'* is the radius of the dividing surface, x^'* and x^'^ are the radiuses of the surface boundaries. 

Assuming that the surface width is small we may expand x around x^. Then the corresponding excesses take the 
following form 

£y [0] {x^) = W}{x') + ^i{(ly}{x') + . . . 

x^ 

<tL[l]{x') ^^{t}{x') 

where '^{q} is the planar excess operator, given by Eq. (|VIII.2p . while '^i{q}{x^) is the first moment of the planar 
excess, given by 

px'^'' - x' 

Viii{q}{x')= axxq^'^ix-^x') (VIII.5) 



We see, that if the radius of the droplet is large, all the excesses may be approximated by the planar interface 
excesses. As the radius decreases, the contributions proportional to 1/x* become more and more significant. Particu- 
larly for small droplets, the resistance has an important curvature contribution. Contrary, the interfacial conductance 
of the droplet does not depend on its size at all. As the spherical droplet has a two-dimensional isotropy, it is also 
follows that the excess conductance is still a scalar for a droplet of an arbitrary size. 



IX. CONCLUSIONS 



In this paper we have build a general approach for the description of the transport phenomena in the interfacial 
region. We have shown how the continuous description can be linked to the excess properties of the surface. This 
makes it possible to consider the surface as an autonomous phase not only in equilibrium but also in non-equilibrium. 

We have derived integral relations for the interfacial resistances for the transport perpendicular to the surface an 
the interfacial conductances for the transport parallel to the surface. They are the coefficients in the linear force-flux 
relations for the heat and mass transfer across and along the interface. The analysis was done for curved surfaces 
which allows us to address different applications where the evaporation and condensation happens. In particular, we 
considered planar interfaces, which are valuable for distillation processes, and spherical interfaces, which are important 
in nucleation processes. 

Building the link between the continuous and discontinuous description allows one to see the important factors which 
affect the transport phenomena most. Among them we can distinguish the local resistivity/conductivity profiles, the 
enthalpy of vaporization and the curvature of the surface. 

The excess of the local resistivity to transport normal to the surface gives the resistance of the whole surface. This 
is similar to the resistance of a finite bulk layer, which may be thought as a series of resistors. Similarly, the excess of 
the local conductivity for transport along the surface gives the conductance of the whole surface. This is also similar 
to the conductance of a finite bulk layer, which may be thought of a set of parallel conductors. If there is a peak in 
the local profile, it immediately affects the resistance of the conductance of the whole surface. It is therefore crucial 
to have the information about the whole profile of a local transport coefficient, not only its bulk values. 

Both the interfacial resistances and the interfacial conductances depend on the variation of the enthalpy across 
the interface. The overall transport coefficients depend on the equilibrium enthalpies which vary a lot through the 



19 



interface. One can see from the above formulae, that the dependence on the enthalpy of evaporation (the difference 
between the enthalpies of the liquid and gas phases) is crucial not only for the diagonal diffusion coefficient, but 
also for the off-diagonal coefficients. This is an important result since cross coefficients are usually neglected in the 
description of the interfacial phenomena. 

The curvature of the surface is important for such processes as nucleation. As the radius of the droplet becomes 
smaller, the higher moments of the resistivities contribute to the overall interfacial resistances. 
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Appendix A: Excess in curvilinear coordinates 



Consider a scalar density 4> being the divergence of a vector function: (j) = V-q(r). Its excess 

(£y[0](a;^r„) = e;y[V•q](x^r||) = -i- / dxi f)i [)2 f}3 (V-q)'^^(r; x^) 
If a density cj) is given by a parallel divergence (/> = Vy -qn (r), then its excess is 



(A.l) 



€y[0](a;^r„) = €y[V||•q,|](x^r||) 



1 



Furthermore, we use the standard formula for the divergence of a vectorial function in curvilinear coordinates 

1 



(A.2) 



( d d d 

f)2 ()3 \axi dx2 0x3 



and the parallel divergence of a vectorial function 

f)i f)2 f)3 \dx2 dx3 
According to the definition pil.4p 

(V-q)'=^(r;a;^) = V-q(r) - {V -(fir)} Q{x' ~ Xi) ~ {V-q^(r)} 8(2:1 - 2;'') 
= V-(q'=^)(r;a;'') + q^'(r)-Ve(a;'' -xi) + q^(r)-Ve(a;i -x") 

where q'^^(r;a;'') is defined similarly to Eq. (jVI.4[) and Eq. pil.4[) . 

Consider V-(q'^^), the first term in the second line of Eq. (jA.5|) . Using Eq. (jA.3[) one can show that 



(A.3) 



(A.4) 



(A.5) 



da:if)if)2t)3V-(q^^)(r;a:^)= f)2[)3 9! 



X"' 



dxi 



(A.' 



^(f)r^3C) + ^(t).^2C) 



The first term in Eq. (|A.6p vanishes since, according to Eq. (jlll.ip . q'^{x^-^) — q'^{x^''^) = 0. It is easy to verify, that 
the second term in Eq. (jA.6p is equal to (}| [}| (By [ Vy -qn ]. 

Consider the last two terms in the second line of Eq. (jA.sp . Using the standard formula for the gradient of a scalar 
function in curvilinear coordinates 



1 dvu 1 dm 1 dvu 

= IT'S — + ir^ — 12 + 7—^ — 13 

f)l dxi [}2 0X2 f)3 0X3 



(A.7) 
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one can show that for Heaviside step function Q 

^ I da;l[)l()2f)3q''(r)•Ve(a;l-a:^) = q''(a:^r„)•il = g^(x^r||) (A.8) 

Eq. (jA.ip then takes the foUowing form 

£y [ V • q ] (x^ r „ ) = (x^ r II ) - (x^ r II ) + £y [ V„ • q,, ] {x' , r„ ) (A.9) 

We note, that Eq. (jA.9|) is the exact equation for the general form of the excess function. No assumptions about the 
curvature magnitude, hke in Eq. (jVI.SP or Eq. (|VI.4p . or the nature of the vectorial function q were made here. 

Appendix B: Excess of a parallel flux for the thin surface 

Consider a scalar function w{x±, ry) such that it is independent on the normal coordinate x±. Furthermore, assume 
that corresponding extrapolated functions va^ and obey the same property. It follows therefore that a function vu 
and the extrapolated function tu'' have the following form: 

tn(a;_L,r||) = ro(r||) 

(B.l) 

ro''(a;_L,r||) = ro^(r||) 

According to Eq. (jlll.ip . ■n7(a;'[''*, ry ) = vj^ {x'];'' , r which leads to 

n7(r||) = tn''(r||) (B.2) 

Consider the excess (By [JirVyZZj], where the function zu satisfies the above relations. Using the definition (|A.2p . 
one can show that 



ev[J,rV||^] = ^y^^^^ dx^ 



0X2 J \ 8x3 



(B.3) 



It follows from Eq. pOj) that 

Using the definition (jVI.Sp . this can be written as 

€v[J||-V||ru] = €y[J||]-V||n7 (B.5) 

Appendix C: Local resistivities 

We need to relate the resistivities r'^ from Eq. pi.l9p to the resistivities r' from Eq. (|VII.6I) . This is done by 
comparing the coefficients at the same fluxes in these equations. To do this we need to translate the set of fluxes used 
in Eq. (IVII.6|) . {J^', Ji, ■ ■ • , Jn-i}, to the set of fluxes used in Eq. pi.l9|) . {Jg, Jji, ■ • • , •/{„}• This is done with the 
help of the relation 



fe=i 

n 
k=l 



(C.l) 



Substituting J' and Ji^^ into the first of Eq. (jVIL6p we obtain 



. n—l n—1 n—1 



k=l k=l 
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Comparing it with the first of Eq. (jII.19[) we obtain 



e ' 
T = T 
' qq ' qq 

n-1 



-r'qqhi-^rljk^k+r'q^, i ^ l,n-l 

Tl-l 

'"''qq^n ~ ^ ^ qk Cfc 



(C.3) 



fc=i 

which are the first 3 equations of Eq. (|VII.7p . 

In order to obtain the remaining relations we consider the second of Eq. pi.lQp , which gives 

j=i j=i i=i j=i 



(C.4) 



i=l 



Furthermore we use Eq. pi.Sp . In case of the transport in the direction only perpendicular to the surface v can be 
taken away and Eq. pi.Sp becomes 

EC.(v|-^.vi)^0 (C.5) 

z=l ^ 

Together with the second of Eq. (|VII.6|) it gives 

n — n 

-Ee^v-f =-E^»^»^4 

^ n-l (C.6) 

- ^ = -/ij„ V_L - + rj^ Jq + ^ r J, 

i=l 

Substituting V(l/T) from Eq. (jC.2p and and Ji from Eq. (jC.l|) we obtain the left hand size of Eq. (jC.6|) expressed 
in terms of the fluxes Je and J^j and the resistivities r' . Comparing the result with Eq. (IC.4P we obtain the following 
equations sets 

(C.7a) 

^j<2 ~ '""9 ~ ~^qq^jn + ''jgj J — l,fJ~l 



ELi '^L-C/c = + el! ^qk ^k ~ r ^) ELi 6 



J* 



(C.7b) 



ELi ^L^fe = + el! <fc ELi ^k hk 



ir'hn + J2Ur' £,k)h 



(C.7c) 



■ el! '^jk ^k - r^gh„, j = l,n-l 
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solving which we obtain the relations between the remaining resistivities 



n-l 

'^I'J = -^qqf^J - ^kq 6 + r-g, j = 1,71-1 

k=l 

fc=l 
n-l 

k=l 

n— 1 n— 1 n— 1 

fe=l 1=1 k=l 
n-l 

= rqqhjhn + S.k{rlqhn + r'^Jij) - rj^/l„ (C.8) 

n— In— 1 n— 1 

k=l 1=1 k=l 
n-l 

fc=i 

n— 1 n— 1 n-l 

fe=l (=1 fe=l 

n-l n — 1 n — 1 

fc=l fe=l i = l 



As one can confirm the symmetry of the r '-matrix leads to the symmetry of the r^-matrix and vice versa. We 
therefore do not give the expressions for rig and r|„ in Eq. dVILTll . 



Appendix D: Local conductivities 



We need to relate the conductivities t"^ from Eq. pi.l8|) to the conductivities I' from Eq. (jVII.Sp . One can do that 
by comparing the coefficients at the same forces in those equations, like it was done for resistivities for the transport 
normal to the surface (see Appendix [C]). We will employ a different procedure for conductivities, however. We will 
compare the coefficients at the same terms in the expression for the entropy production. Both procedures use the same 
method, so they are equivalent. Namely they use the fact that if expressions ^ aiXi — ^ PiXi for linearly independent 
set {xi\, then ai = Pi. Note, that even though the size of the conductivity (resistivity) matrix is different for the two 
sets we are going to relate, all the fluxes and forces in the corresponding set are linearly independent. Thus, both n 
forces {V(l/r), -V(Ati„/r) + /ii„V(l/T), • • • , -V(/z„-i,„/T) + /i„_i,„V(l/r)} are linearly independent, as weh as 
n + 1 forces {V(l/T), — V(/xi/T), • • • , — V(^„/T)} are linearly independent (the similar statement is true for fluxes 
in the case of normal transport through the interface). 

We use the equivalent expressions (IVII.3P and pi. 161) for the local entropy production in the parallel direction. 

We first substitute Eq. (IVII.5[) into Eq. (|VII.3I) in order to obtain the expression for the entropy production in terms 
of thermodynamic forces only, V(l/T) and — V(/i,;„/r) -|-/ii„V(l/r) for i = 1, n — 1. We then substitute /L<j„ = Jii—Jln 
and write this entropy production in terms of the the forces V(l/T) and —V^fim/T) for i = l,n. After some algebra 
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we obtain the following expression for the entropy production 



E 

k=l 



n— 1 71 — 1 n — 1 

1=1 i=l j=l 



j=l 

~n— 1 n— In— 1 



i=i j=i 



E E ( ^" V ' " T 



Lfc=i 



fc=l i=l 



V II „ 



•E 

k=l 



-V 



Vi-l 



T 



V ^ 



EE^i. 

i=i j=i 



(D.l) 



We then substitute Eq. (|IL18[) into pi.l6[) and obtain the expression for the entropy production in terms of the 
thermodynamic forces V(l/T) and — V(^i„/T) for i = \,n. As a result we get 



T 



fe=i 



V- 



^qk ^kq 



EE 

1=1 3 = 1 



T 



(D.2) 



Comparing the coefficients at corresponding terms in Eq. (|D.1[) and Eq. (|D.2I) (in square brackets) we get 



^qq ~ ^qq + (^ig + ^gi) + hinhjn ^ij 

i=l 2=1 j = l 



^q/c + ^feq — ^qk + ^kq + ^'^hin{£lk + ^ki)j k—l,n 1 



^qn ' ^nq 



-I- 



n-1 



n—1 n—1 



T.(Kk + ^k,)~T.T.^^r.{eL + ^L) 



k=l k=l i=l 

iji iji, i,j l,r7, — 1 

n-l 



i=l 
n—1 n—1 

EE^;. 

i=i j=i 



(D.3) 



Without restrictions in generality we assume that the matrix of the coefficients is symmetric, so Eq. (ID.3P can be 
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written as 

n-l 



t'qq — £qq + ^ ] {tj^q + ^qi) + ^ ] ^ ' hinhjn l^j 

i—1 i—1 j—1 



n-l 



i=l 

n — 1 n— 1 n— 1 



fc=l fc=l 1=1 (D.4) 



i,i = l,n-l 

n-l 



J2 J«' 

n-l 



i=l 
n-l n-l 

i=l 3 = 1 



which leads to Eq. (|VIL9[) together with Eq. (|VII.10p . 
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